Multidomain Local Fourier Method for PDEs in Complex Geometries
ثبت نشده
چکیده
but with higher frequencies and amplitudes of the oscillating components), only a low accuracy can be achieved after an aaordable number of PCG iterations. The reason is that a strip-wise domain decomposition cannot handle large changes in non-constant coeecients in the lengthwise direction (see explanation at the end of section (5.1)). More general domain decompositions in the form of rectangular boxes will be able to treat eeciently any complicated geometries. where Laplacian r 2 r;; is written in polar coordinates. This equation is solved repeatedly at each time step by the generalized MDLF algorithm as described in previous sections. The strip-wise decomposition Eq.(5.4) of the transformed (rectangular) domain ! is carried out. For typical used in computations, 1 that justiies the implementation of the local matching procedure for the Helmholtz type equation (5.9), was discussed in section 2. As an example we consider the annular region plotted in Fig.7. The inner and outer boundaries are given, respectively, by The source function in Eq.(5.5) is f(x; y) = e x + e y and the boundary values are (x; y) = e x + e y ; (x; y) 2 @. As t ! 1, the solution u(x; y; t) approaches the steady state regime u 1 (x; y) = e x + e y ; (x; y) 2. The maximum relative error in the steady solution is given in Table 3. It decays exponentially fast as the resolution in the azymuthal direction increases.
منابع مشابه
A New Technique of Reduce Differential Transform Method to Solve Local Fractional PDEs in Mathematical Physics
In this manuscript, we investigate solutions of the partial differential equations (PDEs) arising inmathematical physics with local fractional derivative operators (LFDOs). To get approximate solutionsof these equations, we utilize the reduce differential transform method (RDTM) which is basedupon the LFDOs. Illustrative examples are given to show the accuracy and reliable results. Theobtained ...
متن کاملA Robust Reconstruction for Unstructured WENO Schemes
The weighted essentially non-oscillatory (WENO) schemes are a popular class of high order numerical methods for hyperbolic partial differential equations (PDEs). While WENO schemes on structured meshes are quite mature, the development of finite volume WENO schemes on unstructured meshes is more difficult. A major difficulty is how to design a robust WENO reconstruction procedure to deal with d...
متن کاملComparisons between pseudospectral and radial basis function derivative approximations
Fourier-based pseudospectral (PS) methods have been used since the 1970s for obtaining spectrally accurate solutions to PDEs in periodic geometries. Radial basis functions (RBFs) were introduced about the same time for interpolation on scattered nodes in irregular geometries. As was later recognized, they can also be used for accurate numerical solution of PDEs. Although the main strength of RB...
متن کاملA Hybrid Fourier-Chebyshev Method for Partial Differential Equations
We propose a pseudospectral hybrid algorithm to approximate the solution of partial differential equations (PDEs) with non-periodic boundary conditions. Most of the approximations are computed using Fourier expansions that can be efficiently obtained by fast Fourier transforms. To avoid the Gibbs phenomenon, super-Gaussian window functions are used in physical space. Near the boundaries, we use...
متن کاملImplementation of D3Q19 Lattice Boltzmann Method with a Curved Wall Boundary Condition for Simulation of Practical Flow Problems
In this paper, implementation of an extended form of a no-slip wall boundary condition is presented for the three-dimensional (3-D) lattice Boltzmann method (LBM) for solving the incompressible fluid flows with complex geometries. The boundary condition is based on the off-lattice scheme with a polynomial interpolation which is used to reconstruct the curved or irregular wall boundary on the ne...
متن کامل